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ABSTRACT 

We reexamine stochastic Fermi acceleration (STFA) in the low energy (Newto- 
nian) regime in the context of solar flares. The particle energization rate depends 
a dispersive term and a coherent gain term. The energy dependence of pitch angle 
scattering is important for determining the electron energy spectrum. For scat- 
tering by whistler wave turbulence, STFA produces a quasi-thermal spectrum. 
A second well-constrained scattering mechanism is needed for STFA to match 
the observed 10 — lOOkeV non-thermal spectrum. We suggest that STFA most 
plausibly acts as phase one of a two phase particle acceleration engine in impul- 
sive flares: STFA can match the thermal spectrum below lOkev, and possibly 
the power law spectrum between 10 and lOOkeV, given the proper pitch angle 
scattering. However, a second phase, such as shock acceleration at loop tops, 
is likely required to match the spectrum above the observed knee at lOOkeV. 
Understanding this knee, if it survives further observations, is tricky. 



Subject headings: acceleration of particles. Sun: flares. Sun: X-Rays, gamma-rays 



1. Introduction 



Fermi acceleration was flrst proposed as a mechanism for cosmic ray acceleration (Fermi 
1949, 1954). In the original model, compressive perturbations in the Galactic magnetic fleld, 
associated with molecular clouds, reflect charged particles. If these clouds converge, the 
particles gain energy over time. If they diverge, the particles lose energy. Later, it was 
realized (Bell 1978; Axford, Leer, & Skadron 1978; Krymskii 1977; Blandford & Eichler 
1987), that shock fronts are another site of Fermi Acceleration. In the shock acceleration 
model, charged particles stream into magnetic perturbations in the post-shock region, reflect, 
and are scattered back across the shock by pre-shock Alfven waves. Repeated reflections 
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steadily accelerate particles to a power law distribution. This has been studied extensively 
(e.g. Jones & Ellison 1991). If the shock thickness is determined by the ion gyroradius then 
ions are picked up out of the thermal population but electrons must be injected at energies 
at or above the thermal energy by a factor of the ratio of proton to electron mass, {mp/rrig), 
to incur power law acceleration. The injection process is a critical outstanding problem for 
many applications of shock acceleration. Shock Fermi acceleration is commonly referred 
to as first order Fermi Acceleration because the sign of the energy gain after each cycle is 
positive and dE/dt oc Vc/v, where E and v are the electron energy and speed, and Vc the 
velocity of the magnetic compression. 

Fermi acceleration can also take place in a fully turbulent plasma. There, the mir- 
roring sites are turbulent perturbations, typically fast mode magnetohydrodynamic waves, 
randomly distributed throughout the plasma (e.g. Achterberg (1984)). Electrons encounter 
these perturbations such that there is a stochastic distribution of energy gaining and energy 
losing reflections. As is demonstrated below, there is a net dissipation of turbulent energy 
into high energy electrons. Because the energy gaining and energy losing reflections are equal 
to first order, STFA is often referred to as second order Fermi acceleration and proceeds more 
slowly than the first order process. For STFA, dE/dt oc {ydvY. 

STFA has been considered extensively as the acceleration mechanism in impulsive solar 
fiares, (e.g. LaRosa ct al. (1996)). Observations of these fiares show hard X-ray emission 
with a downward breaking power law spectrum extending from lOkeV to at least O.SMeV 
with the break energy narrowly distributed around Ei^.^. = lOOkcV and a thermal distribution 
at energies below lOkeV (Dulk et.al. 1992; Krucker, & Lin 2002). The time structure of the 
emission shows distinct spikes of duration Is and typical energy lO^^erg (Aschwanden, et.al. 
1995). Non-thermal emission occurs principally at the footpoints of the soft X-ray loop, and 
to a lesser extent at a loop-top hard X-ray source (Tsuneta 1996; Masuda, et.al. 1996). Brown 
(1971) has shown that the emission at a dense target is consistent with Bremsstrahlung 
radiation by electrons accelerated to a power law energy distribution at some height above 
the target; in impulsive fiares the acceleration site can be associated with the loop-top 
region, while the thick target is associated with the footpoints. We demonstrate that STFA 
is possibly responsible for the acceleration of electrons below lOkeV, while first order Fermi 
processes at the loop-top fast shock may produce the highest energy electrons. In this picture 
STFA also provides power law distributed electrons in the range 10keV< E < lOOkeV, thus 
satisfying the shock injection criterion and producing the observed spectral break at lOOkeV. 
It is shown that in order to produce this spectrum the pitch angle scattering must obey the 
restriction that the scattering distance, Ap, is inversely proportional to the energy of the 
scattered electron in the 10 — lOOkcV range. Matching the knee is difficult, requiring either 
a cutoff in the secondary pitch scattering at lOOkeV, or the sudden onset of yet another 
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pitch scattering agent with a much shorter wavelength. While both are possible, these 
requirements provide a serious challenge to STFA models of electron acceleration. 

It is important to note that the downward break is not observed in all impulsive solar 
flares. Indeed, Dulk et.al. (1992) observed a spread in break energies and spectral indexes. 
Some of their flares exhibited no discernible break, or even some upward breaks (ankles). In 
this paper we address the fundamental process of STFA, and show how it can accommodate 
the downward breaking spectra observed in a subset of impulsive flares. The general form 
of the electron spectra we model is illustrated in figure 1. We assume a thick target Bremm- 
strahlung emission model where the electron spectral index 6 and the photon spectral index 
7 are related hj 6 = 7 + I (Brown 1971; Tandberg-Hanssen & Emslie 1988). The data shown 
are mean values taken from the flares studied by Dulk et.al. (1992): Ei„. = lOOkeV, S — A 
below Efyr, and 5 = 5.25 above Ei,r. 

Furthermore, not all flares are observed to be dominated by electron acceleration. A 
recent flare observed by Hurford et.al. (2003) clearly shows regions of emission which are 
dominated by X-ray emission from electrons as well as regions which are dominated by 
gyrosynchrotron emission from MeV/nucleon ions. Proton and ion emission appears to be 
associated primarily associated with larger flare loops. Miller and Roberts (1995); Miller, 
Emslie, and Brown (2004) propose that this can be explained by a two stage process for 
ion acceleration. First, ions are accelerated via gyroresonance to speeds of roughly va by 
Alfven waves, and subsequently are accelerated preferentially over electrons by magnetosonic 
waves. They argue convincingly for the gyroresonant acceleration by Alfven waves on the 
basis of relative ion abundances. However, it is unclear that the second stage acceleration by 
magnetosonic modes must be resonant. It appears that the second stage is consistent with 
STFA. In any event, proton acceleration in long flare loops presents an interesting problem for 
acceleration models, in that suepr-Alfvenic protons must be preferentially accelerated over 
electrons, but electron acceleration still must be dominant in shorter loops. In this work, 
we presume that the loops are sufficiently short that protons remain sub-Alfvenic. Longer 
loops, and the effects of proton acceleration on the shaping of the observed Bremsstrahlung 
X-ray spectra from high energy electrons, require further study. 

STFA is found to depend on two competing effects which we refer to as the steady 
and diffusive acceleration rates. The steady rate represents the net acceleration of electrons 
due to the slight advantage of head-on or energy gaining refiections over catch-up or energy 
losing refiections. The diffusive term represents the spreading of the electron distribution 
function as a result of the stochastic nature of the refiections. Longair (1994) treats the 
two effects together using the Fokker-Planck equation. Likewise, Park and Petrosian (1995) 
discuss general solutions of a simplified Fokker-Planck equation for STFA. In a follow-up 



-4- 



work, Park, Petrosian, and Schwartz (1997) apply their solution to solar flares. While they 
produce spectra consistent with observations, they do not discuss the physics of electron 
escape. Furthermore, they focus mainly on the regime above lOOkeV. Some past treatments 
of STFA in impulsive flares focused exclusively on the diffusive term. LaRosa et al. (1996) 
derives the diffusive acceleration rate using simple physical arguments. Chandran (2003) 
derives the diffusive term using both phenomenological arguments and quasi-linear theory. 
The latter also includes Coloumb losses as a small correction; this is of note since the Coloumb 
loss term is mathematically similar to a negative steady acceleration term. Herein we derive 
the steady term and compare it to the diffusive term, finding the steady acceleration to 
be dominant in the non-relativistic regime for impulsive flares. As will be seen, our steady 
term differs slightly from that of some previous calculations such as that in Longair (1994) 
in its dependence on the turbulent magnetic fluctuation strength. This arises because we 
average only over the pitch angle phase space for which STFA operates in magnetic mirroring, 
whereas Longair (1994) averaged over all pitch angles in considering a more generic form of 
Fermi acceleration. A similar averaging over all pitch angles is performed in Skilling (1975); 
Webb (1983). 

We first review the basic Fermi process using a test particle approach. We then derive 
an expression for the mean acceleration of electrons in a turbulent plasma via STFA, and 
compare this to the diffusive STFA derived by LaRosa et al. (1996). Finally, we discuss the 
trapping of electrons in the turbulent accelerating region and show that at non-relativistic 
energies, the electron spectrum depends strongly on the energy dependence of pitch angle 
scattering. For scattering by whistler wave turbulence, the emerging spectrum is quasi- 
thermal. In order to produce the 10 — lOOkeV power law in solar flares, we show that there 
must be an additional source of pitch angle scattering which has a length scale inversely 
dependent on electron energy; this mechanism is tightly constrained. The existence of such 
a pitch angle scattering mechanism is yet to be determined. This brings to the fore the most 
pressing difficulty with STFA models of electron acceleration; the pitch angle scattering 
requirements are stringent and might not be possible to meet. 



2. The Fermi Acceleration Process 

Consider a particle of charge q traveling with gyroradius tq in a magnetic field of 
strength B (Fermi 1949, 1954; Spitzcr 1956). The charge follows a helical path, orbiting the 
field line while also moving parallel (or anti-parallel) to the field line. Taking the condition 
for circular motion and the Lorentz force 
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where q is the electron charge, and applying conservation laws for angular momentum and 
kinetic energy {rGV±_ and v"^ constant) yields 

qL s\v? 9 , , 

where L is the angular momentum of the charge, E the kinetic energy, and I'sin^ = v±_ 
relates the total velocity to the component perpendicular to the field line. The pitch angle, 
9, is the angle between the field line and the velocity vector. As the charge enters a region 
of increasing such as a magnetic compression, the pitch angle evolves according to 

sin^ 9i sin^ 9 , , 

B + 5B ^ B ^ ' 

where 5B is the increase in field strength and 9i is the pitch angle at field strength B + 5B. 
When sin^i = 1, the charge cannot penetrate further into the compression, and reflects. 
This process, known as magnetic mirroring, is commonly used to confine laboratory plasmas 
(Dendy 1990). It follows immediately that mirroring will not occur at a given compression 
unless the initial pitch angle satisfies 

n B , . 

Fermi (1949, 1954) showed that moving magnetic mirrors, in particular molecular clouds, 
can accelerate charges. In the cloud's frame of reference (primed), mirroring results in only 
a change in the sign of v'^^ , the component of the initial velocity of the charge parallel to the 
field line in the compression's rest frame. Let us work for the moment in the limit where 
the compression speed and the particle speed are both <S c. Transforming to the lab frame, 
5v\\ — ±2i>c, where Vc is the drift velocity of the cloud. The positive (negative) sign is for 
head-on (catch-up) reflections between the charge and cloud. Catch-up reflections are defined 
as those where the components of the compression and charge velocities parallel to the field 
line have the same sign. Head on reflections arc those where the parallel components have 
opposite signs. The net change in energy from a reflection is given by 

5E± = -{vj - v^,) = -{26viivocos{9) + (5^)^) = {2m){±VcVii + v^), (5) 



where Vf and Vq are the speeds after and before reflection, m is the mass of the charge, 
and v\\ — Vocos{9). Head-on reflections result in a positive energy change, while catch-up 
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reflections can result in a negative energy change when v\\ > Vc- In both types of reflection, 
there is the positive term proportional to v'^. 

One can repeat this derivation using Lorentz transformations in place of the Galilean 
transformations to generalize this result to particles of any energy scattered by compressions 
which are still restricted to non-relativistic velocities. For a derivation see Longair (1994); 
we simply cite the result: 

(6) 

where c is the speed of light in vacuum. Over time, charges trapped between converging 
magnetic compressions are subject to only head-on reflections, and are accelerated to higher 
energies. 

Notice that the change in momentum of a Fermi accelerated electron is solely in the 
component parallel to the mean field. This corresponds to an increase in the electron's 
pitch angle. As an extreme example, consider an electron of initial energy Eq with pitch 
angle in the mean field B approaching 7r/2. Upon doubling the electron's energy via Fermi 
acceleration, the pitch angle in the mean field is reduced to 7r/4. Clearly, acceleration to 
high energy must be accompanied by some additional scattering agent which isotropizes 
electron pitch angles on a short time scale, otherwise Fermi acceleration shuts off after small 
accelerations as pitch angles evolve out of the range given by (4). The well known problem 
of pitch angle scattering remains largely unsolved (e.g. Achterberg (1981); Melrose (1974); 
LaRosa et al. (1996)). 

Fermi acceleration is distinct from the transit time damping (TTD) treated, for exam- 
ple, ill Miller, Larosa, & Moore (1996). TTD is the magnetic analog of Landau Damping. 
In TTD, electrons (or ions) which are near gyroresonance with waves of wavenumbwer k 
are pushed towards the resonance by field gradients in the wave which alter the parallel 
component of the velocity. Gradients in the electron velocity spectrum result in a net damp- 
ing or enhancement of the waves. In the presence of a spectrum of waves, an electron can 
drift from resonance at k to k ± 5k and so forth, eventually reaching high energies. Fermi 
Acceleration, however, is a non-resonant interaction. Electrons will mirror at a compression 
regardless of energy provided that the pitch angle is sufficiently large. TTD is often referred 
to as resonant Fermi Acceleration because the two processes rely on similar physics. Table 
1 lists the relevent length and time scales for STFA in impulsive flares. 



SE±R = 2E 
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Length scale 



Description 



Time scale 



Description 




Turbulent outer scale 
Turbulent eddy scale 
Parallel eddy scale 
Perpendicular eddy scale 
Effective STFA scale 



Ted 



V 



Steady STFA acceleration time 
Diffusive STFA acceleration time 
Turbulent eddy time 
Pitch angle scattering time 



MFP for pitch angle scattering 
Ap for scattering by whistlers 

Constrained A„ for 10 - lOOkoV 



Table 1: Table of length and time scales relevant to STFA in impulsive flares. 



3. 



Stochastic Fermi Acceleration 



We now consider the behavior of charges in a turbulent magnetic plasma where magne- 
tosonic modes provide the sites of magnetic mirroring. This scenario differs from first order 
Fermi acceleration by shocks in two ways. 1) Consecutive mirroring events are not coher- 
ent, but rather stochastically distributed between head-on and catch-up. 2) The turbulent 
cascade governs the acceleration efficiency; the system picks out a scale where acceleration 
competes with the cascade. In the solar corona plasma, Vc, the velocity of the magnetic 
compressions, is the phase speed of the magnetosonic modes, which is roughly the Alfven 
speed for the fast mode and the sound speed {cg) for the slow mode. Typically, thermal 
electrons in the corona are super-Alfvenic and non-relativistic, va ^ vq <^ c, and (3 ^ 0.05; 
we will solve the STFA problem in this regime. 



Recall that the energy gain from a typical reflection is given in Eq.(5) to be 5E± — 
2m{±v\\Vc + v'^). We define three parameters: R, the total rate of refiections; the rate 
of head-on refiections; and the rate of catch-up refiections. The relation R — R+ + R^ 
is automatically satisfied by this definition as all reflections must be of either the head-on 
or catch-up type. This allows us to write the approximate acceleration rate as the sum of a 
coherent term and an incoherent term: 



3.1. Determination of the Steady Acceleration Rate 
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where the subscript S is used to distinguish our derived acceleration rate from that of LaRosa 
et al. (1996). The first term, proportional to (i?+ — i?_) represents the mean acceleration 
due to the offset in the rates of the two types of reflection. The second term, proportional to 
(i?+ + i?_) = R, the total reflection rate, represents the coherent term. This expression gives 
a full description of the mean acceleration of charges by the non-relativistic STFA process. 
To evaluate it in a particular plasma requires the determination of the head-on and catch-up 
reflection rates i?+ and 

In order to obtain and consider the path a charge takes to encounter a mirror. 
Note that there is a distinction between an encounter and a mirroring because of the pitch 
angle condition for reflection. We take the fraction of encounters which reflect to be F and 
assume that this fraction is the same for both head-on and catch-up encounters: F — — 
F_. This assumption is often taken in the regime where v Vc (LaRosa et al. 1996). While 
this assumption is not strictly true, the effects of relaxing it are negligible. In Appendix B 
we repeat this calculation without assuming = F_ . There is a well defined mean distance 
between encounters. At as well as a relative velocity between the particle and the compression 
v± = v\\ ± Vc- The small difference in the head-on and catch-up speeds is responsible for 
the offset in rates. For each type of encounter, the mean separation is 2Xt and the rate of 
reflections of any sort is specifled by the general relation R = Fv^eiK with v^ei the relative 
velocity between the particle and compression, yielding 



V\\ -\- Vc 

R+ = F- " 



2A 



T 

'c 



ii_ = F!l^ (8) 

V\\ 

R^R+ + R-^F-^. 

At 



The offset in rates is thus 



i?+-i?_ = F-^, (9) 
At 

and we can rewrite the average acceleration rate as 

/ dE\ V\\v'^ V\\vJ. vnvt , , 

— = 2mF-p + 2mF-p = 4mF^, (10) 
\ dt J g At At At 

with the associated acceleration time scale ts — E / We thus see that the acceleration 

due to the offset in head-on and catch-up reflection rates and the acceleration from the 
coherent term in Vc are equal. 
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To simplify the derivation of the electron spectrum, we recast (10) as 

dE\ , dE dM 

— = ^mvlR = -r^—r. (11) 
dt ) g " dM dt ^ ' 

where we used (8) and define M as the total number of refiections experienced by an electron, 
dM/dt = R, and 

dE/dM = Amvl. (12) 

The quantity dE/dM, the mean acceleration of an electron per reflection, will be of 
particular use when examining electron escape in section 4. 



3.2. Comparison to the Diffusive Acceleration Rate 

A different approach was taken by LaRosa et al. (1996). They set = i?_ in the 
V Vc regime and studied the diffusion of particles through energy space via random walk. 
The starting point in their calculation of the electron acceleration was the timescale for the 
e-folding of a charged particle's energy in the turbulent plasma 



1 1 fdE 



td E\dt J ^ NSt' 



(13) 



where is the number of mirrorings required to double the particle's energy and 6t is the 
time between encounters. They set R^ — /?_ = 0, and also dropped the last term in in 
Eq 6. From (7) and (10) it is clear that if one of their two assumptions is valid the other 
must also apply, and {dE/dt)s = in that limit. Under these assumptions 5E+ and 5E- 
are equal in magnitude, and from the standard solution of an evenly weighted random walk 
N — {E/5Ey. The acceleration rate is then 

where Vo is the total initial speed and the subscript D is used to denote LaRosa et al's 
diffusive acceleration rate. 

To complete the calculation of the acceleration rates, we must obtain F and average 
over pitch angles. The minimum accessible pitch angle for reflection is related to F by 

F = cos(^^i„)= (— j , (15) 
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where we have applied the reflection condition from (4). In the regime where the compression 
ratio 5B/B « 1, which will apply to the plasma of interest, and taking the assumption 
that pitch angles arc isotropic gives (cos(6')) = cos(^^j„)/2 and (cos^(^)) = cos^(^^j„)/4. We 
now write the averaged acceleration rates as 



What do these two acceleration rates represent? {dE/dt)s is the steady growth of the 
mean kinetic energy due to the drain of turbulence by the combined effects the slightly non- 
zero {R+ — it!_) and the coherent v\ term: a shift of the mean electron energy to higher 
energy. On the other hand, (dE/dt)]:) represents the diffusion of energies away from the 
mean via random walk: a spreading of the distribution. The relative importance of the two 
is fixed by their ratio 



The combined result of the action of both processes on an initially narrow Gaussian energy 
distribution is shown in Fig. 2 where we have chosen ( = 65. and assume that electrons 
do not escape. Thus we can examine the evolution of electron energy spectra solely due 
to the inffuence of the two acceleration rates. As ( is increased, the steady (mean growth) 
term becomes increasingly dominant over the diffusive (distribution widening) term. To 
understand STFA in a particular plasma, both acceleration rates must be calculated. In the 
event that the diffusive rate is very small compared to the steady growth rate, it can be 
ignored. The steady growth rate is always faster than the diffusive growth rate for a > 1. 
As will be shown later, ( ^ 100 in flare plasmas, and the diffusive term is neghgible. 

It is very important to note that our result differs from the standard for Fermi Accelera- 
tion, in which both the diffusive and steady terms depend on the same power oi Sb/B. This 
difference arises as a result of the averaging over pitch angles. To correctly obtain (AE), 
one must only average over those encounters which result in a reflection. For traditional 
STFA, the range of pitch angles which reflect is ultimately determined by the turbulent mag- 
netic field strength. One factor of cos(^) in the expression to be averaged results in one 
factor of (dB/By^"^ in the acceleration rate. In other treatments, such as that of Longair 




(16) 




(17) 
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(1994); Webb (1983); Skilling (1975), the acceleration mechanism is assumed to act at all 
pitch angles. In this case, the averaging over cos{6) while maintaining the asumption of pitch 
angle isotropy, yeilds a numerical value with no 6B dependence. In this regime, the diffusive 
and steady acceleration terms have the same relative strength at all levels of turbulence. 

3.3. Specification of the Turbulent Cascade 

This leaves SB/B as the remaining parameter to be determined. It is related to the tur- 
bulent length scale At through the cascade law. Magnet ohydro dynamic (MHD) turbulence 
proceeds by the shredding of like sized eddies and subsequent formation of smaller eddies; en- 
ergy input into eddies on a large (outer) length scale, Lt, cascades rapidly to smaller length 
scales on the eddy turnover time, T^d, and finally dissipates at A^, the dissipation scale. If the 
cascade obeys Komolgorov's steady state assumption then the energy flow through all length 
scales is constant, and independent of the scale. This results in an inertial range between the 
scales Lt and A^ where the turbulent energy density has a power law dependence on eddy 
size. The draining of turbulence at eddy size A^ is usually determined by a micro-physical 
process, such as resistivity. When STFA is active, the turbulence can instead be drained 
by pumping energy into electrons. This sets another condition for STFA to proceed in a 
plasma: there must be some Xsf which is greater than the resistive length scale at which 
the STFA timescale is shorter than the cascade time, otherwise the turbulence will drain at 
the resistive scale before STFA can produce an appreciable electron acceleration. 

There are three major MHD wavemodes: Alfven waves, and the fast and slow magne- 
tosonic waves. Alfven waves are purely transverse, and thus do not compress the magnetic 
fleld; they cannot participate in STFA. Both the fast and slow modes are compressive, and 
are in principle capable of Fermi acceleration. It has been argued that in low f3 plasmas such 
as the solar corona, the slow mode is rapidly dissipated via Landau damping (Achterberg 
1981). However, more recent studies of MHD turbulence indicate that the cascade time for 
GS turbulence is significantly shorter than the electron damping time, and slow mode damp- 
ing by electrons can be ignored in turbulent flare plasmas (Lithwick and Goldreich 2001). A 
key difference in the two analyses is that Lithwick and Goldreich (2001) treats MHD turbu- 
lence as inherently anisotropic, whereas Achterberg (1981) assumes isotropy. Furthermore, 
Maron (private communication) has shown in simulations which neglect damping that the 
slow mode may be driven with much higher total energy content than the fast mode at 
low p. A deflnitive resolution of the issue is beyond the scope of this paper. However, we 
should point out that up to this point, the calculation is independent of the choice of wave 
mode. There is one signiflcant difference between the two, however: slow modes propagate 



at roughly the sound speed {vc ^ c^) while fast modes propagate at roughly the Alfven speed 
(vc ^ va)- It will be shown that due to the influence of pitch angle scattering, STFA is likely 
dominated by the fast mode. 

MHD turbulence is in general anisotropic; the direction of any large scale mean magnetic 
field defines a preferred axis. Also, even if the turbulence is isotropic on large scales, smaller 
scales may see the larger scale turbulent structures as an effective mean field. Goldreich 
and Sridhar (1997) (hereafter, GS) modified the the Komolgorov assumption for the cascade 
of slow and Alvfen modes of MHD turbulence under the condition that the turbulence is 
anisotropic with scale Ay along the field line and Aj_ perpendicular to the field line. The two 
directions are found to obey different cascade laws, with Ay cascading more weakly than A^. 
The parallel direction is of more importance to STFA, as it represents the distance along the 
field fine between refiection sites. The GS power law for the parallel scale is (Goldreich and 
Sridhar 1997; Lithwick and Goldreich 2001) 



where B is the mean magnetic field strength, and 5B is the turbulent field strength at parallel 
length scale Ay. 

The exact power law of MHD turbulence remains the subject of some debate, so for 
now we assume a general power law of form 



where a > 1 is an arbitrary index. Using the turbulent power spectrum and substituting for 
F from (15) we can rewrite the acceleration rates as 




(19) 




1/a 



(20) 




(21) 



For a typical turbulent cascade, where A^ < Lt and a > 0, ts < td- 



3.4. XsF and the Role of Pitch Angle Scattering 

There remains the final step of determining the particular dissipation scale A 5^7 at which 
the energy drain takes place. In order for STFA to overcome the cascade of MHD turbulence, 
it must drain energy at a rate equal to the input rate at the outer scale. If turbulent com- 
pressions are shredded and cascade faster than electrons can draw out energy via reflections, 
then the cascade continues to smaller length scales. At smaller scales, STFA is more rapid. 
STFA becomes competitive with the cascade at a scale determined by {dE/dt)s — dET/dt. 
{dE/dt)s is acceleration rate of electrons and dEx/dt is the cascade rate of turbulent energy, 
nnipV^/ Lt- The STFA scale Xsf is different for acceleration by the two compressive MHD 
modes. 

For slow mode turbulence, Vc — Cs, and the balance is 

2nme 2 f \ ^^^""^'^ uTJip 3 
Solving for Xt/Lt and associating this particular At with Xsf gives 

^SF 

In solar flares and a GS turbulent cascade {a — 2), vq — 1.2 x lO^cm/s, = 3 x lO^cm/s 
and va — X lO^cm/s (LaRosa et al. 1996), this gives Xsf/Lt ^ 10~^. We have taken 
the initial electron velocity to be the mean velocity of the thermal background plasma. The 
cascade will proceed down the inertial range to this length scale where STFA then acts as the 
micro-physical damping agent, rapidly draining the energy from turbulence into particles. 

In the case of the fast mode, where Vc — va, the rate balance is 

2nme 2 f \ nnip 3 

and the STFA length scale is then given by 

Lip 

For solar flare conditions, and a GS cascade (a = 2), Xsf/Lt — 10~^. We have tacitly 
assumed that the length scale for pitch angle isotropization is roughly equal to At- If it is 



(22) 



1 rrip 

2 rUe civ 



(23) 



(24) 



InipVA 
2me V 



(25) 
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not, the acceelration rate is retarded significantly, and STFA can be shut off. To understand 
this we must further explore the role of pitch angle scattering. 

As discussed above, pitch angle scattering is necessary during acceleration to maintain a 
population of electrons which satisfy the pitch angle condition for refiection. The strength of 
the pitch angle scattering strongly regulates the rate of acceleration. We consider three cases: 
^ XsF, Ap <S XsF, and Ap ^ A^^ where eddy \p is the typical distance over which pitch 
angles are isotropized. In the first case, Ap ^ ^sf, electrons refiect a few times and quickly 
leave the pitch angle range in which they can refiect. They then must stream a distance 
of order Ap before they can scatter again. Thus the rate of refiections and the acceleration 
rate are both decreased by a factor of A^^/Ap. Since the acceleration rate and cascade rate 
are not in balance, the cascade continues down to smaller scales Ay < Xsf- The nominal 
acceleration rate (eq 22) is proportional to Xrp^^^, while the retardation factor is proportional 
to At. The combined effect is a net acceleration rate which is proportional to X^^; smaller 
scale turbulence is actually less efficient as an accelerator. As a result, STFA never turns on 
in this regime. In the second case, Ap ^ Xsf, the pitch angle scattering is far more rapid 
than acceleration. Since pitch scattering can take an electron through /i = 0, very strongly 
pitch scattered electrons traverse the plasma by random walking in steps of length Ap, again 
reducing the rate of refiection, this time by a factor of (Ap/A^i?)^. Unlike the previous case, 
this is not a problem for STFA; the retarding factor tends towards unity as the cascade 
continues to scales At < Xsf- The net acceleration rate is proportional to A^^^^, and STFA 
turns on as the cascade proceeds to a sufficiently small scale. In case three, where Ap ^ A5, 
pitch angle scattering and refiections proceed at the same rate. Thus, electrons are capable 
of streaming freely from compression to compression, while they maintain a nearly isotropic 
pitch angle distribution. This is the simplest pitch angle scattering regime for STFA. 

The identity of the accelerating wavemode is now easy to determine. In section 4.1 we 
show that whistler wave turbulence is a plausible source of pitch angle scattering, at least 
for the lower energy quasi-thermal component of the spectrum. At 3keV, Xwh/Lx is roughly 
10^"^. This places slow mode turbulence {Xsf/Lt ^ 10~^) well in the first regime. Slow 
modes do not participate in STFA in these fiares. Fast modes, however, have Xsf/Lt 10^^ 
and therefore are in the nearly ideal range for acceleration. Furthermore, both Xsf and Ap 
grow linearly with electron energy, so as electrons undergo STFA by fast modes in the 
presence of whistler wave turbulence, they remain in the same pitch angle scattering regime 
throughout. 
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4. The post-acceleration spectrum 

The simplest case of STFA is the steady state, where we assume that electrons are 
injected into a turbulent region at energy Eq at a rate equal to that at which accelerated 
electrons escape. The turbulent energy supply is continuously replenished at a large scale. 
We are concerned with the energy spectrum, N{E), of electrons escaping the region. Note 
that this is in general different from the spectrum of the electrons within the turbulent 
region. We define Nt{E) to be the total number of electrons reaching energy at least E 
before escaping, such that 

N{E) (X -dNt{E)/dE. (26) 

Initially, we consider the case of strongly relativistic electrons; a full derivation of this regime 
is presented in the Appendix. 

To appreciate the calculational differences between the non-relativistic regime of inter- 
est to solar flares and the more commonly studied relativistic regime, we begin with the 
latter. Following the approach used by Bell (1978) for Shock Fermi acceleration, and writing 
dNt/dM — —pesc^t, where Pesc is the mean probability of an electron escaping from the 
acceleration region, gives 



dNt dNtdM ^dM ^1 , , 

lE = -dMlE ^ -P-^'lE ^ -P-^'^^ (2^) 
where on the right hand side we have for the moment taken the strongly relativistic limit: 
dE/dM — aE and assumed Pesc to be constant. This treatment of the highly relativistic 
limit follows Fermi (1949). One can solve for Nt{E) by separating variables, integrating both 
sides and inverting the logarithms, resulting in the familiar power law (see e.g. Fermi (1949); 
Longair (1994); Jones (1994)) 

Pesc 

Nt{E) = N(^^j " , (28) 
where is the total number of electrons. From Eq. (26), one obtains 

N{E) = No i^—j , (29) 

where N^dE is the number density of escaped electrons with E — Eq. Notice that the 
logarithmic integrals in both Nf and E are vital to producing the power law. 

For STFA by fast mode waves, the computation is more comphcated because Pesc is 
energy dependent. We solve for a general Pesc in the non-relativistic regime, leaving the 
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specification of the trapping for later discussion. In the non-relativistic regime, the accelera- 
tion rate is not proportional to the kinetic energy as it is in the strongly relativistic regime. 
Instead, one has dE/dM — Amv\. We thus write, using (11), 

dNt _ NtPesc s 
— 2 ' 



dE Amv\ 



which can be rewritten as 



^ = ^^dE- (31) 

We have assumed that dN^jdM = Pesc^t- This is reasonable as long as the electrons can be 
treated as statistically independent and coUisionless. In this case, it is reasonable to assume 
that Pesc carries no inherent dependence on Nt and the escape rate is simly given by the 
product of the number of electrons in the volume and the mean escape probability. Taking 
Pesc — Po{E/Eq)~^ allows us to solve for a particularly interesting N{E). We see immediately 
that 

= -^J^dE, (32) 

and N{ET) is again a power law energy distribution: 



= ^0 1^;^ J , (33) 

where b = poEo/Amv\ . However, in any other case, STFA does not produce a simple 
power law. The importance of the trapping mechanism is now clear; the combined energy 
dependence of the acceleration and escape must be E~^ to produce a power law spectrum. 
Such a spectrum rehes on the coincidental logarithmic integrals over both N and E. 



4.1. Calculation of Pesc 

Let us now calculate p^sc for non-relativistic electrons within the turbulent volume, and 
consequently the energy spectrum of electrons. For simplicity, let us take the turbulent 
region to be rectangular, with the long axis, z, parallel to the direction of the bulk flow, 
with z — and z — Lp fixed to the downstream and upstream boundaries of the region 
respectively. L^^ is taken to be the extent of the region of turbulent flow, which is presumed 
to be the entire distance between the reconnection sheet and the top of the soft X-ray loop. 
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This distance is typically of size 10^° cm for solar flares (Tsuneta 1996). The largest eddy size 
in the turbulence, Lt is set by the width of the outflow, typically lO^cm. Thus the turbulent 
volume consists of a number of cells, each of which flows downward from the reconnection 
point towards the loop-top. An electron escapes the acceleration region only when it reaches 
the X-ray loop at the base of the turbulent region. These individual cells may be associaetd 
with single bursts or fragments of X-ray emission, and thus are responsible for the temporal 
structure of impulsive flares. In order to escape the region with energy E{M), an electron 
must stream from its location in the region at some height z to the boundary at z = 
after the Mth reflection without further reflection. We will assume that the electrons are 
contained in the region in the x — y plane by gyration around large scale field lines. To further 
simplify the problem, we shall assume that the electron density remains uniform throughout 
the turbulent region. We also neglect the bulk flow speed, = 8 x 10^ cm s^^ (Tsuneta 
1996) since the legth of the downflow region is roughly 10^'^cm. This gives a flow time from 
the reconnection region to the loop-top of 100s. The acceleration process is fixed to the much 
shorter Is time scale by the temporal size of the observed energy release fragments and the 
MHD eddy turnover time. Thus, bulk flow into the flare loop is not hkely to be a dominant 
process in cutting off the acceleration. 

Take the mean ^-component of the distance streamed between reflections to be A^; 
carries an energy dependence inherited from the energy dependence of the pitch scattering. 
The probability of escaping at z = after the Mth reflection from a point at height z is 
given by 

Vesc{z) = ^e-^/'% (34) 

and the mean escape probability of electrons distributed uniformly across the length of the 
region is 

Peso = ^ rPesc{z)dz = ^(1 - g-^^^/^^), (35) 
J^F Jo J-iF 

where we have taken A^ = Ap/2 from the isotropy in pitch angles, with \ the pitch angle 
scattering length scale. To obtain the spectrum of solar flare electrons requires specification 
of the pitch angle scattering. 

Both Miller, Larosa, & Moore (1996) and LaRosa et al. (1996) assume strong scattering, 
and suggest that the scattering agent above IkeV is resonant interaction with lower hybrid 
(LH) turbulence, or circularly polarized electromagnetic waves, such as whistler waves. Below 
IkeV, Coulomb interactions are thought to be sufficiently strong to isotropize the electrons. 
It should be noted that Melrose (1974) sets the threshold for the whistler mode resonance 



at 25 keV for flare plasmas, while Miller & Steinacker (1992) argue that the resonances 
extend down to IkeV. We assume the latter. In a recent study, Luo, et.al. (2003) considered 
whether LH wave turbulence is the primary mode of electron acceleration in solar flares. 
They concluded that the pitch angle scattering is too inefficient to maintain isotropy. Thus, 
we assume that LH wave turbulence cannot supply sufficient pitch angle scattering to sustain 
STFA either. Whistler waves are more promising. 

Melrose (f974) associates the frequency of pitch angle scattering v with the pitch angle 
diffusion coefficient in the quasi-linear equation. Thus, f/z/ is the chracteristic time scale for 
effective pitch angle isotropization of the electron distribution. It should be noted that, in 
general, some small anisotropy is likely to remain in the distribution, and that this anisotropy 
could be responsible for the generation of the whistler waves. However, the source of these 
waves is still uncertain. Prom Melrose (1974), we have that 



y= ' ^ \ , (36) 

where 7e ^ 1 is the Lorentz factor of the electron, is the electron gyrofrequency, cUp is the 
plasma oscillation frequency, and 

nip 

is the energy density of the turbulence at the resonant wavelength. This allows us to rewrite 
(36) as 



i/ = 5 X lO^Sioonio^/^ 



^'•^ , (38) 



where we have used v\ = B'^/ATTUp, = 1.8 x lO^-Bioo, oOp = 5.7 x lO^n^g^, and the 
dimensionless parameters i?ioo = -B/IOOG and nio = ng/lO^^cm"^. 

We must compare u to the growth time for pitch angle anisotropy due to STFA. Bearing 
in mind that for STFA, dv±/dt = 0, 



dE dv\\ 



and the pitch angle evolves according to 



d{cose) 1 sin'^e dE 



dt 



m cos Q dt 



(40) 
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Substituting in from equation (21)for dE/dt and assuming a GS cascade (a = 2), gives 

d{cose) 2 (\t\~'''^v\ 1 

- -ItT -tATt) ^^^^ 

For impulsive flares, Eq — Q2keV , B^qq — 2, and n^Q — 1. Taking = \sf results in 
i/ = 1 X 10^ s"^ and 

.2 n^o\ „_i 



= 2 X 10^ (^-^J (42) 

where vsf is evaluated at the threshhold pitch angle for reflection. To maintain pitch angle 
isotropy, scattering by whistler waves must occur on a time scale shorter than pitch angle 
evolution by STFA. Thus, as long as Usf < isotropy can be maintained. This condition 
is met for all E > Eq. Whistler modes, if present, are capable of providing sufficient pitch 
angle scattering to maintain isotropy. 

In addition to maintaining pitch angle isotropy, the scattering mechanism must also 
operate at a length scale which traps electrons in the volume; X^h <^ Lt = lO^cm must 
be satisfied, or else electrons rapidly leave the acceleration region and STFA shuts off. We 
obtain the pitch scattering length scale for whistler waves, Xu,h, 



1/2 

A^ft = ^ = 2x10^1^—) cm. (43) 



The required condition is satisfied for energies below lOOkev. 



4.2. The Electron Spectrum and Constraints on the Secondary Pitch Angle 

Scattering 

In order to obtain the spectrum of the escaped electrons, we can now substitute the 
functional form for Pesc from (35) into (31) 



dN^ ^ Peso ^ _ e-WAp)^ 

Nt Amv\ Lf Amvjs^ 



^dE^^n- e-^^/^") ^dE. (44) 



Next, by choosing Ap = A^/j, we use \p/ Lp = AE^^"^, with E in units of keV, and rearranging, 
obtain 
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N{E) 



dNt 



1 



(45) 



(IE 



4mv\ 



The resulting spectrum, N{E), is plotted in figure 4. This is consistent with the thermal 
component to the flare spectrum observed using RHESSI (Krucker, & Lin 2002). 

In addition to the thermal component, RHESSI observations show a clear power law re- 
gion extending from roughly lOkcV up to at least 50kcV, above which the data arc uncertain, 
but consistent with a continuing power law. Previous observations using the ISEE3/ICE in- 
strument also show a power law throughout the range of the instrument, 25 — 300keV; the 
spectrum typically breaks downward at lOOkeV (Dulk et.al. 1992). More recent observa- 
tions with RHESSI could push the low energy threshold for the power law as high as 35keV 
(Holman et. al. 2003). The spectral index below the break is ^ 3, while above the break 
it is ^ 4. Elsewhere (Blackman (1997); Selkowitz and Blackman (2004) in preparation), we 
discuss first order acceleration at the loop-top fast shock. Fast shocks are well known to 
accelerate super-thermal particles to power law energy spectra, even in the non-rclativistic 
limit (Bell 1978). However, in order to be accelerated, electrons must satisfy the requirement 
that E ^ {mp/mejVg = lOkeV in solar flare plasmas, where Vg = lO^cm s~^ is the inflow 
speed of the plasma at the shock (Blackman & Field 1994). This places the injection energy 
at rougly lOOkeV. The correspondence of the shock injection energy and the observed break 
energy is noteworthy. A possible mechanism to reproduce the observations is for STFA to 
produce a power law spectrum in the 10 — lOOkeV regime which then satisfies the injection 
criterion for loop-top fast shocks. Since STFA by magnetosonic turbulence in the presence of 
whistler wave turbulence pitch scattering is insufficient to produce the power law component. 
However, it is possible that a second pitch scattering agent exists which produces the power 
law in the 10 — lOOkeV range. We examine the constraints imposed on this scattering. 

To produce a power law spectrum, non-relativistic STFA requires Pesc oc E~^. From 
(35), we see that this is true only if \p/Lp = T/E, where 2E /T ^ 1 and the exponential 
term is small. Here F is a constant parameter which fixes the strength of the pitch scattering. 
While the physics of the pitch angle scattering is not well understood, this constrains the 
scattering mechanisms available to STFA. We define Ac = LtV/E — 2 x 10^(£'o/£')cm to 
be the pitch scattering length scale of the constrained mechanism, and the electron spectum 
is given by 



r is constrained by the observed X-ray spectral index of 7 = 3. It is a standard prediction of 
flare models (Brown 1971; Stepanov and Tsap 2002; Kiphnger, et.al. 1984) that the spectrum 




(46) 
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of electrons accelerated above the loop-top is steeper than the spectrum of the thick target 
Bremstrahlung X-rays emitted at the footpoints in solar flares. We assume that to match 
the RHESSI X-ray data requires an index ^ 4 for the electrons, or 

r = 12mvl = 0.072A;ey = 0.17£;o. (47) 

The exponential term in p^sc is indeed small as 2E/r — 270 at £^ = lOkeV. The electron spec- 
tral index could conceivably be as high as 6, in which case F = 0.28£'o, and the exponential 
can still be safely neglected. 

It is insufficient to merely produce the proper power law. The scattering agent must 
also be able to reproduce the transition from thermal to power law spectrum at the correct 
energy, Ec- We recall X^h to be the length scale of pitch angle scattering associated with 
whistler wave turbulence. In impulsive flares, Xp = Xwh below Ec- Above E^ Xp = Xc- To 
obtain E^ one sets X^h — Xc- Ec — 23keV, which is consistent with the observed threshold 
of - lOkeV. 

There is one more important constraint imposed by the observations: the knee at 
lOOkeV. Dulk et.al. (1992) demonstrate a distinct downward break in the power law spec- 
trum at roughly E^r = lOOkeV. The break energy varies somewhat from flare to flare, but is 
consistently observed in all of the impulsive flares in their sample. Unlike downward breaks, 
upward breaks are easily explained by the meshing of two acceleration mechanisms, as the 
shallow component which dominates above the break emerges naturally from beneath the 
steeper power law which dominates below. For upward breaks, E^r is the naturally occurring 
crossover point. The matching problem is much more difficult for knees in the absence of 
significant cooling on timescales of interest. Since the steep component is above the break 
energy and the shallow component below it, both must be truncated at the break energy. 
If either one extended beyond the break, then that one would overrun the other, and there 
would be no break at all. The most natural solution for a knee is a single acceleration mech- 
anism which undergoes some transition at the break energy. One such example is the knee 
found by Bell (1978b) in the spectra of shock accelerated electrons at roughly IGeV. This 
knee results from the transition from the non-relativistic to the relativistic regime. There 
is no apparent natural transition for STFA of electrons at lOOkeV. However, there is a well 
defined low energy cutoff for a power law spectrum at lOOkeV, the shock injection energy. 
The shock injection threshold is not only at the right energy, but is also a variable cutoff, 
depending on the ion temperature and local magnetic field strength, consistent with the 
variability in the observed Ef,r. 

Bell (1978b) has shown that shock acceleration does not change the spectrum of electrons 
if the pre-shock spectrum is shallower than the post-shock spectrum which obtains from a 
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steep pre-shock spectrum. Shock Fermi acceleration cannot steepen a power law spectrum; 
it can only make it shallower. This is another difficulty for knee matching. The STFA 
spectrum must have a sharp cutoff at Ei,r in order to match the knee. This may not be an 
impossible condition to meet, especially as Ei,j. may be greater that the injection energy, not 
precisely equal to it. 

One natural cutoff occurs when = A^; acceleration will shut off when the cascade 
reaches the resistive scale. For slow modes in impulsive flares, A^ = lO^cm (LaRosa et 
al. 1996; Chandran 2003; Lithwick and Goldreich 2001), which for F = 0.073keV gives a 
cutoff energy of 7 x 10''^ kcV, which is both too high and far outside of the non-relativistic 
regime. A possible solution is that the constrained pitch scattering has a maximum resonant 
threshold at E^r- Another possible solution is that yet another very strong pitch scattering 
mechanism has a threshold energy of Ef^ and a length scale < A^. Both of these solutions 
are presently ad hoc. This underscores the need for a more thorough understanding of pitch 
angle scattering in astrophysical plasmas. It also illustrates the limitations of STFA as an 
acceleration mechanism in solar flares; if STFA alone were to account for the spectrum from 
10 — 100 keV, the tight constraints on the pitch angle scattering mechanism that we have 
identified are required. 



5. Summciry and Discussion 

STFA in the non-relativistic limit behaves differently from highly relativistic STFA. 
At the core of these differences is the energy dependence of the electron velocity at low 
energies. Thus, unlike the relativistic case, both the rate of reflections and the probability 
of escaping the acceleration region at an energy E vary. Using a test particle approach, we 
have examined this behavior and derived the spectrum of post-acceleration electrons in a 
plasma under impulsive solar flare conditions. 

For traditional STFA, where there is a minimum pitch angle constraint which determines 
whether an individual encounter results in reflection, it is seen that the steady acceleration 
rate can dominate over the diffusive acceleration rate. This arises from the averaging over 
pitch angles to evaluate (AE'). Some previous treatments of the generahzed Fermi accel- 
eration problem do not have such reflection conditions, and thus do not retain factors of 
the turbulent field strength, SB/B, when averaging. In those treatments, such as Longair 
(1994); SkiUing (1975); Webb (1983), the steady and diffusive terms typically are seen to be 
of the same order. For some processes this is appropriate, however non-resonant STFA is not 
one of them. Thus, the phase space conditions for scattering by the acceleration mechanism 
can play a very significant role, even in cases where pitch angle isotropy is maintained. 
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The nature of the pitch angle scattering turns out to be the dominant factor in deter- 
mining electron escape, and therefore the shape of the spectrum. We find that whistler wave 
turbulence, which is well studied in solar flares (Melrose 1974; Miller & Steinacker 1992), is 
an excellent source of pitch angle scattering which allows STFA to produce a quasi-thermal 
electron distribution that peaks at E ^ 5keV. This matches the lowest energy portion of 
the observed X-ray emission very well. However, to produce the power law spectrum ob- 
served in the range ^ 10 — lOOkeV by STFA requires at least a second scattering mechanism. 
Matching the spectral index and the transition energy from quasi-thermal to power law 
spectrum requires an undetermined scattering mechanism which satisfies Xc/Lp — V/E 
with r = 0.073A;ey, and naturally becomes the dominant pitch angle scatterer at roughly 
20keV. 

If the constrained pitch angle scattering mechanism is discovered, it implies that the 
acceleration of electrons in solar fiares is at least a two stage process. The first stage, STFA 
in the downflow region, produces both the quasi-thermal spectrum below ^ IQkeV and the 
lower half of the power law spectrum up to lOOkeV. To produce the highest energy electrons, 
as well as the spectral break dX E = lOOkeV requires a second acceleration mechanism at 
the top of the soft X-ray loop. We are further exploring the possibility that first order 
acceleration at a weak fast shock, formed as the downfiow impacts the top of the closed flare 
loop, is responsible for electron acceleration to the highest observed energies. Acceleration 
at fast shocks is known to have an injection energy of roughly lOOkeV, and varies with 
temperature. This coincides with the break in the power law spectrum at lOOkeV, and is 
consistent with the variabihty observed by Dulk et.al. (1992) in E'^^. 

Recently, Chandran (2003) concluded, using quasi-linear theory, that STFA for slow 
modes is not viable in the 10-lOOkeV regime. While we also find slow modes to be ineffective, 
differences between our paper and Chandran (2003) must be kept in mind. Chandran (2003) 
assumed that dp/dt oc p for STFA. While this is true in the strongly relativistic limit for 
STFA, we do not assume that this is true in the lower energy regime (see eq. 22). Second, 
unlike Chandran (2003) we do not assume herein that Pgsc has to be energy independent. 
These two assumptions play a significant role in shaping electron spectra. 

Another concern which can be raised about the effectiveness of STFA as the electron 
acceleration engine in impulsive fiares is the total energetics of the process. Since STFA, as 
developed above, only is efficient in a short length scale regime where one also has 5B/B ^ 1 
it might seem that only a small fraction of the released flare energy is available for electron 
acceleration. This is not the case. The total energy contained in single turbulent cell is given 
by {l/2)mpv\nL^ ^ lO^^erg, where n — 10^°cm~^ is the electron number density in the flare 
plasma. The energy in a single turbulent cell is similar to the energy contained in one X-ray 
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emission fragment. Although only a fraction of the energy in one turbulent cell is ever at Xp 
at one time, it does all cascade down to Xp over an eddy turnover time. Thus, while SB/B 
is always small, the energy throughput can still be high enough to accelerate the electrons. 
The similarity in total energy between a single turbulent cell and an individual impulsive 
X-ray fragment strongly suggests that the two are related. 

Miller, et.al. (1997) estimates that as much as ~ 94% of the magnetic energy in a flare is 
available in the turbulence, which is sufficient to produce the high energy electrons inferred 
from the observed X-rays, but raises concerns about the efficiency of STFA, particularly in 
competition with other sources of dissipation. While we have not fully studied other dissipa- 
tion mechanisms which might compete with STFA for this energy, three significant ones can 
be ruled out: proton acceleration by STFA, Landau damping, and resistive dissipation of the 
turbulence. The latter two have already been discussed. Proton acceleration is a significant 
concern since Fermi acceleration of protons and heavy ions was in fact the very problem 
Fermi intended to solve. Thermal protons in coronal flare plasmas are sub-Alfvenic and thus 
cannot meet the condition for mirroring (LaRosa et al. 1996; Blackman 1999). However, 
Miller and Roberts (1995) argues convincingly that gyroresont interaction of protons and 
heavy ions with Alfven waves can accelerate them to velocities above on a relatively short 
time scale. Within their model, the ions then are accelerated by compressive magnetosonic 
waves at the expense of electron acceleration. Recent RHESSI observations (Hurford et.al. 
2003) indicate that the emission signatures of ions and electrons are spatially separated, with 
the ion emission associated with longer loops. Miller, Emsfie, and Brown (2004) concluded 
that these observations are consistent with the gyroresonance model of ion acceleration; as 
the loops grow longer, protons are more likely to reach super- Alfvenic speeds, and thus can 
be accelerated by the magnetosonic waves. This second phase of acceleration need not be 
gyroresonant. While it appears promising, further study is required to determine if STFA 
models can accomodate these results. 

The strong dependence of the post-acceleration electron spectrum on the pitch angle 
scattering agent is both a positive and negative feature. It leaves STFA considerable flexibil- 
ity in matching various characteristics of solar flare X-rays which fall outside of the simple 
scenario studied in this paper. For example, Lin et.al. (1981) flrst observed a superhot com- 
ponent in a solar flare, which has since been supported by RHESSI observations Krucker, 
et.al. (2003). This thermal, or nearly thermal, spectral component is seen at energies of up 
to 35keV. Within our STFA framework, the superhot emission can easily be explained by an 
enhancement of pitch angle scattering at lower energies, either by increased whistler wave 
turbulence, or some other scattering agent. While this flexibility naturally allows for the 
wide range of flare characteristics observed, it does not yet deflnitively solve the flare accel- 
eration problem. Instead, it shifts the focus exclusively to a well constrained, but largely 
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unspecified, array of pitch angle scattering mechanisms. This is the single greatest obstacle 
to STFA models of acceleration. 

In short, STFA can naturally account for the thermal spectrum below lOkeV, and some- 
what less naturally for the non-thermal spectrum between lOkeV and lOOkeV. There we have 
shown that Xp must depend inversely on particle energy, in contrast to that of pitch angle 
scattering by whistler waves below ^ lOkev, which is proportional to the particle energy. 
Above lOOkeV, shock acceleration is a natural possibility; the needed injection of super- 
thermal electrons may be provided by STFA operating at energies below Ebr- The knee at 
lOOkeV remains the most difficult spectral feature to accommodate, and we have explained 
the difficult requirements to pitch angle scattering that this demands. 
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Appendix A: The Power Law Spectrum of Relativistic STFA 

Notice that the spectrum we obtain for STFA is different from the power law result 
of Jones (1994). This is a matter of regime; we discussed in the text the acceleration of 
non-relativistic particles in a region of non-relativistic turbulence, here we show that when 
the particles are relativistic, a power law spectrum emerges. 

Recall that (5) for fully relativistic electrons in a region of non-relativistic turbulence is 
given by (6) 



SE±R = 2E 



where E is the total energy, kinetic plus rest, of the electron before reflection. Notice that if 
we the low velocity limit, f ^ c where E = mc^, the expression reduces to (5). The relative 
velocity between the compression and electron for head-on and catch-up type interactions 
are still given by (8) so the steady acceleration rate is given by 



dE\ 



2E 



{R+ - R. 



c c 



A 



(48) 



Alternatively, we can find the mean acceleration per refiection by multiplying equation 48 
by 
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In the highly relativistic hmit, E is just the kinetic energy, and we recover the famihar result 
(Jones 1994) that dE/dt oc E. This proportionality is expected to produce a power law. 
We derive the power law spectrum for STFA of highly relativistic electrons by following the 
approach of Bell (1978) and assume that Pesc is a constant in flare plasmas, independent of 
electron energy. We start by integrating dE/dl to obtain E{1) 



l^^\n{E/E,), (50) 



Av 



2 



-,2 



The probabihty of an electron remaining in the acceleration region for at least I reflections 
is given by 

P(/ + ) = (1-Pesc)'. (51) 

Taking the logarithm and substituting in for 1 from equation 51, gives 



I / E\ f E^ (-pe.cM)-l 



In P{E+) = Hn (1 - p,,,) = - In j In (1 - p,,,) = In j . (52) 

where we used the approximation ln{l — Pesc) = —Pesc for Pesc <^ 1 in obtaining the expres- 
sion to the right of the final equals sign. Differentiating with respect to E, results in 



P{E) = Eo l^-j , (53) 

where P{E)dE is the unnormalized probability of a post-acceleration electron having the 
energy E. In the limit, where Pesc is extremely small, the relativistic STFA spectrum has 
power law index ~ 1. In a plasma where Pesc ~ ^, the power law index can grow larger, and 
the index is very sensitive to Pesc- In the third regime, where Pesc ^ electrons stream out 
of the turbulent volume quickly, do not experience much acceleration, and have a very steep 
power law energy distribution with virtually no very high energy electrons {E ^ Eq). 
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Appendix B: Derivation of Steady Acceleration Rate with F+ ^ F_ 

In section 3.1 we derived the steady acceelration rate for electrons in a low /9 turbulent 
magnetic plasma. This derivation was contingent on the assumption that = F_ = F, 
which is not strictly valid. Blackman (1999) calculates F for Fermi acceleration. By resetting 
the limits of the integral in his eq (12), and renormalizing for the smaller phase space, one 
arrives at 



where cos 0^ is the minimum pitch angle at which an electron will reflect and f ^/f is the ratio 
of the Alfven speed to the electron speed. We rename these quantities A and B respectively; 
both are small quantities. By taking a series expansion of eq (54) and truncating it at second 
order in S, it can be simplified to 




(54) 



A 1±B--B^ 



(55) 



2 



Recall that from (8), 




(56) 



Prom this one easily obtains 




(57) 



and 




(58) 



This gives us all of the ingredients for calculating the steady acceleration from (7) 




The resulting acceleration rate is 
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AB^ 2A + A'^--B- -AB"^ + , 



(59) 
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A 



where we have added the additional subscript h to indicate the distinction from the previously 
calculated rate. The steady acceleration rate found in (10) from the assumption = F_ = 



Note that provided A > AB this is the largest term in (59). Indeed, for coronal flare plasma, 
B ^ 0.1 at electron energy Eq and decreases with increasing energy while A ^ 0.1 as well 
at Eq, but is largely insensitive to eletron energy. At the onset of the power law regime, 
E = lOkeV = 30£'o, and B ^ 0.01; all terms of order B^ or higher can be neglected, as can 
the term in A^. Thus we can safely use the assumption F_|_ = F_ = F in this regime, and 
(11) is reasonable. 
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Fig. 1. — The typical downward broken power law of an impulsive flare, as observed by 
Dulk et.al. (1992). The left panel shows the hard X-ray spectrum, with Ei^ — lOOkeV, and 
spectral indices above and below Ejjj. of 4.25 and 3 respectively. The right hand panel shows 
the electron spectrum in the emission region inferred from the given photon spectrum using 
a thick target Brcmstrahhmg model for the emission (Brown 1971; Tandberg-Hanssen & 
Emshe 1988). Again, Ebr — lOOkev and the spectral indices above and below E^r are 4 and 
5.25. 
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Fig. 2. — The evolution of a sample electron energy distribution with an initially Gaussian 
velocity distribution. The peak of the distribution is evolved from the flare thermal energy, 
0.2keV to the post-STFA mean energy at 16keV. A)The initial distribution function. B) 
The same distribution after being evolved only by the steady process. C) The distribution 
evolved through both the steady and diffusive processes. Note that the relative width of the 
electron energy distribution, IS.E jEm dereases with increasing mean energy Em- 
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Fig. 3. — From (Blackmail 1997). Sketch of a typical impulsive solar flare. Note that the x- 
point reconnection occurs in the filled region at the top of the diagram. Only the downward 
half of the outflow is shown. Reproduced by permission of the AAS. 
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Fig. 4. — The spectrum of non-relativistic STFA under impulsive flare conditions with 
whistler wave turbulence as the only pitch angle scatterer. Eq — 0.3keV. 
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Fig. 5. — The spectrum of non-relativistic STFA under impulsive flare conditions with 
whistler wave turbulence and a second source of pitch angle scattering. The second pitch 
angle scattering source obeys the constraints required to produce a power law. £'0 = 0.3keV. 



